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ABSTRACT 

Context. Twisted magnetic fields should be ubiquitous in the solar corona, particularly in flare-producing active regions where the 
magnetic fields are strongly non-potential. The magnetic energy contained in such twisted fields can be released during solar flares 
and other explosive phenomena. It has been shown recently that reconnection in helical magnetic coronal loops results in plasma 
heating and particle acceleration distributed within a large volume, including the lower coronal and chromospheric sections of the 
loops. Hence, the magnetic reconnection and particle acceleration scenario involving magnetic helicity can be a viable alternative to 
the standard flare model, where particles are accelerated only in a small volume located in the upper corona. 

Aims. The key goal of this study is to investigate the links and observational signatures of plasma heating and particle acceleration in 
kink-unstable twisted coronal loops. 

Methods. We use a combination of MHD simulations and test-particle methods. These simulations describe the development of kink 
instability and magnetic reconnection in twisted coronal loops using resistive compressible MHD, and incorporate atmospheric strat¬ 
ification and large-scale loop curvature. The resulting distributions of hot plasma let us estimate thermal X-ray emission intensities. 
The electric and magnetic fields obtained are used to calculate electron trajectories using the guiding-centre approximation. These 
trajectories combined with the MHD plasma density distributions let us deduce synthetic hard X-ray bremsstrahlung intensities. 
Results. Our simulations emphasise that the geometry of the emission patterns produced by hot plasma in fiaring twisted coronal 
loops can differ from the actual geometry of the underlying magnetic fields. In particular, the twist angles revealed by the emission 
threads (soft X-ray thermal emission; SXR) are consistently lower than the field-line twist present at the onset of the kink-instability. 
Hard X-ray (HXR) emission due to the interaction of energetic electrons with the stratified background are concentrated at the loop 
foot-points in these simulations, even though the electrons are accelerated everywhere within the coronal volume of the loop. The 
maximum of HXR emission consistently precedes that of SXR emission, with the HXR light-curve being approximately proportional 
to the temporal derivative of the SXR light-curve. 
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1. Introduction 

Solar flares are energetic phenomena commonly understood as 
fast releases of magnetic energy stored in the solar corona. The 
coronal plasma heated during such events produces characteris¬ 
tic emission signatures in the soft X-ray range, and the particles 
accelerated in the flaring regions produce non-thermal emission 
tipically in the hard X-ray energy range (see, e.g., the review 
by |Benz||2008} [Fletcher et al.||20lT] ). Twisted magnetic flux- 
ropes can store significant amounts of magnetic free-energy in 
the solar corona, and are susceptible of developing magneto- 
hydrodynamical instabilities which lead to the release of this 
energy. Coronal loops undergoing a kink instability, in partic¬ 
ular, go through an initial exponential growth phase (the linear 
phase of the insta bility), until they start reconnecting with the 
background field ( [Browning et al.|[T008| ). They then relax into 
a state with lower magnetic energy and less twist (hence re¬ 
leasing a fraction of the magnetic free energy stored initially). 
It has been shown recently that magnetic reconnection in heli¬ 
cal (twisted) magnetic structures may result in impulsive plasma 
heating capable of producing the main properties of thermal X- 
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ray emission in solar flares ( [Pinto et ar][2015[ ), and in particle 
acceleration distributed within large volumes of plasma (includ¬ 
ing the lower corona and the chromosphere) leading to non- 
thermal X-ray emission ( [Gordovskyy et al.[[2012[ [2ofil [2014 ). 
The kink instability, therefore, may be a viable alternative to the 
standard flare model, which assumes that particles accelerated to 
high energies in small volumes of plasma near the flaring loop 
apexes are transported downwards along the loop legs, and that 
these interact with the denser chromosphere producing hard X- 
ray emission there and providing hot plasma which would there¬ 
after fill up the loops and emit thermally in soft X-rays. One 
long-standing difficulty the kink instability scenario faces is that 
its triggering requires that the pre-flare loops are strongly twisted 
(see, e.g, discussion in [Bareford et al.|2()T^ , which observations 
of flaring coronal loops rarely seem to indicate (for both confined 
and ejective flares). Notable exceptions exist, however, such as 
reported by [Srivastava et al.| ( [2010[ ), who measured a total twist 
angle of about 127r on a twisted coronal flux-rope before the oc¬ 
currence of a B5.0 class flare (see also simulations of this event 
by [Botha et al.|2012| ). [Pinto et al.[ ( [2015[ ) pointed out that, how¬ 
ever, the morphology of the thermal emission during the initial 
phases of a flare may lead to an important underestimation of the 
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actual twist of the underlying magnetic field lines. Furthermore, 
the kink-instability scenario, along with other scenarios resulting 
in distributed acceleration/re-acceleration ( [Turkmani et al.|200^ 
[Cargill et al.||2006|) may allow one to circumvent t he so-called 
number problem ( |Brown|p~976t [Brown et al.||2009| ), one of the 
main difficulties of the standard flare model. 


In the present study, we focus on deducing observational 
signatures inherent to the triggering of the kink instability in 
twisted coronal loops, combining thermal and non-thermal ef¬ 
fects (plasma heating and particle acceleration) with application 
to confined flares. We will consider here loop models consisting 
of twisted magnetic fiux ropes embedded in a gravitationally- 
stratified background (comprising a hot corona, a transition re¬ 
gion and a fraction of the chromosphere), which represent small 
(or moderate) confined fiares. The origin of such twisted coronal 
loops can be thought of as the result of photospheric source ro¬ 
tation or (see e.g. [Brown et al.|2003][Dalmasse et al.|2014|) or of 
flux-rope emergence through the photosphere (e.g. [Luoni et ak] 
[2011[ [Archontis et al.[|20T^ [Pinto & Brun[[2013] ). But the de¬ 
tailed analysis of the formation of such twisted coronal loops is 
out of the scope of this study; we only discuss here the evolution 
of fiux-ropes already twisted and close to the onset of the kink 
instability. 

We compare light-curves and time-dependent intensity 
maps of thermal (free-free continuum) and non-thermal 
(bremsstrahlung caused by energetic electrons) radiation in the 
keV and deca-keV ranges. We investigate how the spatial distri¬ 
bution of the emitted photon fluxes relate to the dynamical and 
geometrical properties of the simulated loops, how the total con¬ 
tinuum emission spectra evolves in time and on how the prop¬ 
erties of the emission measures respond to the plasma heating 
processes occurring in the magnetic loops. 


2. Model description 

We study the evolution of a twisted coronal loop embedded 
in a magnetised and gravitationally-stratified background at¬ 
mosphere. The coronal loops are twisted magnetic fiux-ropes 
with large-scale curvature (the loops are nearly semi-circular), 
and with both foot-points anchored in the chromosphere. The 
thermo-dynamical evolution of the system is calculated us¬ 
ing fully compressible resistive MHD by means of three- 
dimensional numerical simulations. The MHD computations 
provide the basis for the estimation both of the thermal emis¬ 
sion by the plasma (in the soft X-ray range), and of the hard 
X-ray emission produced by charged test-particles (electrons). 
The methods employed are described in detail in the following 
sections (Sects. |j!T][2.2[ and [Z^ . 


2.1. MHD equations and model parameters 


The evolution of the twisted loops is determined by solving the 
Eqs. W in a 3D (v,y,z) domain taking into account atmo¬ 
spheric stratification, thermal conduction, and radiative losses.fg 
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Fig. 1. Plasma temperature T and numeric density as a function of 
height during the preparatory thermal relaxation phase at two different 
positions above the surface. Solid lines correspond to the initial state. 
Dotted, dashed and dot-dashed lines represent the temporal evolution 
of these quantities during the thermal relaxation phase at equally spaced 
time intervals after the initial state.The two top panels correspond to the 
n and T profiles along a vertical line right above the loop foot-points, 
where the magnetic field is mostly vertical. The two bottom panels show 
the same quantities along a vertical line which crosses the centre of the 
loop, where the magnetic field is mostly horizontal. 


The system is described in terms of magnetic field B, plasma 
velocity v, density p and specific energy e: 
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where p is plasma resistivity and the current density is j = 
V X B/po, and g is the gravitational acceleration, which is con- 
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Fig. 2. Temperature distribution shortly before the onset of the kink 
instability. Upper panel (a) shows the temperature distribution in the 
mid-plane x = 0, while the lower panel (b) show the difference be¬ 
tween temperature distribution shortly before the instability and the ini¬ 
tial temperature distribution. Selected magnetic field lines are shown in 
green at panel (a). Light blue lines at panel (a) denote cross-sections at 
which the temperature and density distributions are shown in Figures]^ 
and|7] Black dashed lines at panel (b) denote AT = 0 levels. 
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stant throughout the whole domain. Plasma pressure and tem¬ 
perature can be expressed, respectively, as p = {y-\)pe and 
r = (y - l)fhlkB€. Here y is the ratio of specific heats of the 
plasma, m is average particle mass (which is half the proton 
mas rrip, corresponding to fully-ionised hydrogen plasma). The 
electric field, required for the test-particle modelling, can be ex¬ 
pressed as 

E = -V X B rjj . (5) 

The Acond. Qjouie. 2visc and Lrad terms in Eq. account for vari¬ 
ations of specific energy due to thermal con^ction. Joule heat¬ 
ing, viscous heating and radiative losses, respectively. Thermal 
conduction is defined as 


S,ond = '^-[{-'<oT^'^B-VT)B] + V- 


-Ko 


7 ^ 5/2 _ 


Bi 


mi n 


vr 


( 6 ) 


where ^min is some very small parameter. Hence, with non¬ 
zero magnetic field {B » ^min), we get Braginskii conductivity 


along the magnetic field, but negligible in perpendicu¬ 
lar direction, while with negligible magnetic field {B ^min) it 
is and isotropic. The following characteristic dimensional 

factors are used: = 4.47 x 10“^ T, po = 2 x 10“^^ kg • m“^ 

and Lo = 10^ m, yielding the characteristic values for Alfven 
velocity vq = 2.82 x 10^ m • s“\ time = 0.36 s, temperature 
To = 3.23 X 10^ K, and electric field Eq = 1.26 x 10^ V • m-^ 
The Joule heating is Q-^ouie = where p is the local resis¬ 

tivity p = T/back + ^crit- The background resistivity 77back is con¬ 
stant throughout the simulation domain and represents classical 
resistivity. Its value is 10“^ ^^pQpQL^B~^ (corresponding to the 
Lundquist number 10^), i.e. this resistivity component does not 
make any noticeable contribution to the system evolution. The 
critical resistivity is defined as 

r]^i = [ 10 "^ ylmpoLlB-\ j >jcnt ^ 

\ 0^ j < Jcrit 


where jcrit = jiP/Po VT /Tq. This component represents anoma¬ 
lous resistivity with the threshold similar to that for ion-acoustic 
instability, when the current drift velocity exceeds local thermal 
velocity. The value of ji constant is chosen so that before the 
kink instability pcrit = 0 everywhere in the domain. The cho¬ 
sen resistivity values have no exact physical justification, as in 
most other MHD models. However, the chosen values appear 
to be realistic, as they provide realistic reconnection times (tens 
of second s). The radiative losses are calculated using the func¬ 
tion from |Klimchuk et ar] ( |2008| ). The shock viscosity and cor¬ 
responding heating rate used in our simulations are described in 
|Arber et aL] ( [200 1| ). The adopted viscosity values are not rep¬ 
resentative of real viscosity in the coronal plasma, and the cor¬ 
responding viscous heating rate is usually lower than the ohmic 
heating rate due to anomalous resistivity. However, high velocity 
gradients during the fast reconnection phase can sometimes re¬ 
sult in substantial viscous heating, comparable to or even higher 
than ohmic heating. Ultimately, the viscous effects result in con¬ 
version of magnetic energy to internal energy of the plasma and 
can be considered as an additional resistivity effective at loca¬ 
tions with high velocity gradients (see Bareford, et al, submitted, 
for more detailed description of this effect). The above equations 
([T]|^ along with Eqs. |5|6|an m. are solved using the lagrangian 
remap code LARE3D ( jArber et al.|200l| ). We consider four mod¬ 
els: three with small coronal loops with different magnetic field 
magnitudes (models C, S and V), and one case with a large coro¬ 
nal loop (model Y). The simulation boxes are cubes with the 
size of 20 Mm in “small” models (C, S and V) and 80 Mm in 
the “large” model Y; see Tablefor details. The numerical grid 
consists of 256 X 256 X 512 elements with uniform step in each 
direction (i.e. the resolution in z-direction is twice as high as the 
resolution in the horizontal directions). 

All the models undergo an initial preparation phase (con¬ 
sisting of a thermal relaxation stage, and a footpoint twisting 
stage). The initial magnetic field geometry is that of a poten¬ 
tial field arcade in a gravitationally-stratified atmosphere. As in 
[Gordovskyy et al.j ( |2014| ), the initial atmosphere is constructed 
of two isothermal regions: a lower region with thickness of 
about 2 Mm representing the chromosphere with temperature 
< lO"^ K and density exponentially decreasing from 5 x 10^^ to ~ 
10^ cm“^, and an overlaying corona with temperature -0.8 MK 
and nearly constant density ~ 10^ cm"^. This idealised initial 
configuration is not in thermal equilibrium, however, due to the 
presence of magnetic field-aligned thermal conduction. We first 
let the system relax until it reaches quasi thermal equilibrium. 
Eigure[2 shows how the vertical profiles of density and tempera- 
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Fig. 3. Magnetic field topology and current distribution at 4 representative instants 
of the simulations (from top to bottom: t = 0, 30, 120, 196 s after the onset of the 
kink instability; model S), illustrating the onset and growth of the kink instability, the 
reconnection episode that follows, and the relaxation of the coronal loop towards a non- 
twisted state. The lines on the left column are magnetic field-lines, with the blue lines 
having seed points very close to fiux-rope’s foot-point centres. The axis are labelled in 
units of Mm. The right column shows iso-surfaces of current density j - jcru (see Eq. 
[^and the discussion which follows). 



Fig. 4. Evolution of the magnetic, internal and ki¬ 
netic energy in model S after the onset of the kink 
instability. 



Model C 

Model S 

Model V 

Model Y 

Box size, Mm 

20 

20 

20 

80 

Bfp, IQ-‘^T 

45 

90 

220 

665 

Blt, 

5 

10 

24 

67 

Rpp, Mm 

0.7 

0.7 

0.7 

2.8 

Rlti Mm 

2.4 

2.4 

2.4 

9.6 

/, Mm 

24.6 

24.6 

24.6 

97.0 


Table 1. Pre-instability parameters of twisted loops in different models. 
Magnetic field values Bfp and Bu are absolute magnetic field values 
at the foot-points and loop-top, respectively; Rpp and Rpr are the loop 
cross-section radii at the foot-points and loop-top, respectively, I is a 
loop length. 


ture evolve during the relaxation process. It can be seen that the 
temperature transition around z = 2 - 4 Mm becomes gradually 


Article number, page 4 of 


14 


smoother. The temperature distribution in the vertical direction 
tends approximately to T (z) ~ (z + const, where the mag¬ 
netic field-lines are vertical.The equilibrium temperature distri¬ 
bution will, however, be considerably different where the field is 
nearly horizontal (as in the middle of the domain, across the cen¬ 
tre of the loop arcades). Heat transfer during the relaxation phase 
results in some evaporation: plasma from nearly vertical field in 
the transition region (z = 2 - 3 Mm) flows upwards, noticeably 
increasing the loop plasma density up to z = 10 Mm. After about 
2000 ^ 0 , temperature and density distributions become nearly sta¬ 
tionary (the resulting atmosphere is thermally relaxed). 

Foot-point twisting applied to the initially untwisted poten¬ 
tial field in order to build the unstable twisted loop configura¬ 
tions. The twisting phase does not represent an actual rotation of 
photospheric foot-points (an effect which is often observed, but 
out of scope of this study). We can use here rotation velocities 
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tribution These circular vortical motions are applied at the 
zones of maximal (the foot-points), and thus are spatially co¬ 
incident with the contours of the vertical field, which are approx- 
imately circular ther e {cf. discussions by |Torok & Kliem|2003| 
[Aulanier et al.|2Q05| ). 

Such foot-point rotation effectively results in continued he- 
licity injection with the total twist increasing by ~ 0.015 tt/s, i.e 
after about 500 s the total twist is of about S/r. At this stage, 
the twisted loop arcade is slightly thicker and taller than ini¬ 
tially (due to the increase in magnetic pressure within the mag¬ 
netic loop) but its large-scale structure stills maintains its quasi¬ 
circular symmetry. The system is close to the onset of the kink 
instability (but still in mechanical equilibrium), and we define 
this state as the initial condition {t = 0) hereafter. The configura¬ 
tion of the magnetic flux-rope, the plasma density and tempera¬ 
ture at this stage are represented in Fig.[^ and the corresponding 
pre-instability parameters are listed in Table HI 

-10 -5 0 5 10 



X 

Fig. 5. Foot-point twisting profile and vertical magnetic field ampli¬ 
tude at the lower boundary (for model S). The grey-scale shows the 
amplitude of the vertical magnetic field B^ (white and black represent¬ 
ing positive and negative polarity) and the blue lines are streamlines of 
the surface motions (see Eq.[^. The axis are labelled in units of Mm. 


considerably higher than those observed in the photosphere in 
order to speed up the simulations, as long as these remain sub¬ 
stantially lower than the local Alfven speed. Under these condi¬ 
tions, the twisting phase should not have noticeable dynamical 
effects on the kink instability and reconnection (see, for more 
detail, |Gordovskyy et al.|2Q14[ Bareford, et al, submitted). The 
foot-point centres are located at [0.0, +6.5 Mm, 0.0] in models 
C, S and V, and at [0.0, ±26 Mm, 0.0] in model Y. Foot-points 
are rotated with the following profile 

1 - tanh \-(r - Rfp)/Kaii\ 

0 J(r) = ojo -^, (8) 

where r is the distance from a foot-point centre, ojq ^ 0.02 s“^ 
is the frequency of rotation, Rfp is the foot-point radius, which 
is 0.7 Mm in models C, S and V, and 2.8 Mm in model Y, 
and L^aii = 0-13 Mm is a scaling factor. Taking model S as 
a reference case, the foot-point rotation period is about 300 s, 
with a corresponding velocity at the edge of the foot-points of 
1.4 X 10"^ m/s. The coronal Alfven speed is close to 7 x 10^ m/s, 
and the Alfven travel time along the whole loop is ~ 35 s. 
Hence, foot-point twisting motions do not generate transients 
which propagate into the corona, and the coronal field has am¬ 
ple time to adjust and remain in near-equilibrium as it is being 
twisted. The other models (V and Y) have stronger magnetic 
fields and higher corresponding Alfven speeds, for which quasi¬ 
steadiness is even better achieved. Chromospheric Alfven speed 
is, of course, much lower, and approaches the footpoint rotation 
speed. However, the chromosphere has a plasma yS » 1, with 
the sound speed remaining considerably higher than the twisting 
speed throughout the whole chromospheric layer (up to one or¬ 
der of magnitude higher for model S). These reasons altogether 
ensure that the system evolves through a sequence of equilibria 
(even though small amplitude dynamical features, such as waves 
can be seen during the driving period). But most importantly, 
the forcing (twisting) time-scale is always higher than the kink 
instability growth-time. 

Figure shows the streamlines of the rotating motions ap¬ 
plied at the surface together with the vertical magnetic field dis- 


2.2. Synthetic soft X-ray emission from a mode! loop 


We calculate the thermal X-ray continuum emission from the 
model’s number density n and temperature T at each point of 
the domain, as in [Pinto et al.| ( |2015] ). The plasma emissivity at a 
given photon energy hv is defined as 

e (hv, T) = eon^T~^^^gff (hv, T) exp j > (9) 


where gff (hv, T) is the Gaunt factor for free-free bremsstrahlung 
emission. The coefficient eq is 6.8 x 10“^^ if th e emissivity is to be 
expressed in units of erg • cm“^ • s“^ • Hz“^ ( Tucker||l975 1. The 
Gaunt factor is approximated by the piece-wise approximation 


gff(hy,T) = 


1 , 

\ hv ) ’ 


hv < kiyT 
hv > k}jT 


( 10 ) 


The photon flux density emitted at the photon energy hv is de¬ 
fined as 


where EM = is the emission measure of a finite volume V 
of plasma of density n and temperature T. The coefficient /q is 
1.07 X 10“"^^ for a photon fiux measured at a distance of 1 AU in 
units of photons • cm“^ • s“^ • keV“^ The total photon fiux over 
a given spectral band is computed by integrating Eq. O over 
the corresponding range of values of hv. We compute the photon 
fiux at different photon energies for each individual grid cell {i.e, 
volume element), each one having well-defined values for the 
emission measure, density and temperature. As the corona is op¬ 
tically thin to X-ray radiation, the total fiux emitted is obtained 
by adding the individual contributions over the whole loop (or 
over a region of interest). 

We estimate the distributions of the differential emission 
measure DEM {T) in our simulations by computing the emission 
measure of the plasma regions whose temperature lies within 
successive temperature intervals at a given time and normalis¬ 
ing it by the temperature interval width AT. This distribution 
then quantifies the amount of plasma at different temperatures, 
independently of the temperature interval considered. That is, 
the differential emission measure are defined here as a function 
of temperature as 

= (12) 
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Fig. 6. Temperature distributions during and after the kink instability. 
Three loop cross-sections are shown (corresponding to the blue lines 
in Fig. 1^: left foot-point (z = 1.9 Mm) (left column), left loop leg 
(z = 4.25 Mm) (middle column) and vertical plain crossing the loop top 
(right column). Times shown are (after onset of the kink instability): 
panel (a) - 0 s , (b) - 81 s, (c) - 196 s, (d) - 377 s. 


where the index k runs through all the plasma elements (grid- 
cells in the simulations) which lie within the temperature inter¬ 
val [T, T + ST]. The values of % and dV^ are the number density 
and the volume of each element. In other words, we first compute 
a temperature histogram with a given temperature bin size AT. 
Then, we verify which grid-cells have a temperature T within 
each of the bins and we sum over all the corresponding individ¬ 
ual EM. Variations in density for plasma at a given temperature 
are implicitly accounted for. Total emission measures EM can be 
obtained by integrating Eq. in respect to temperature, in the 
relevant interval of temperatures. 


Fig. 7. Logarithm of density distributions. Locations and times are the 
same as in Figure 


tion ( [Gordovskyy et al.|2010[|Gordovskyy & Browning|2011| ): 

dr 
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d{yv\\) 

dt 


dp 

dt 
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tn (yv||)^ 
q jk^B 
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dp 

dt 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


coll 


2.3. Test-particle trajectories and synthetic hard X-ray 
bremsstrahlung 


Electron distribution functions are calculated using the GCA 
test-particle code based on first-order guiding-centre approxima- 


Here r, u and vy are the particle gyro-centre position, transver¬ 
sal drift velocity (_L B) and longitudinal velocity (|| B), respec¬ 
tively; Z? = I is the magnetic field direction vector, Ue = 
is the E X B drift component; p is the magnetic moment per 

mass unit p = where Ug is the particle gyration velocity. 
The relativistic coefficients are defined as y = (1 - 

and K = Jl - where the absolute particle velocity is 
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V = + Wg. Coulomb collisions (introduced through 

[^]co// taken into account through averaged de¬ 

celeration and pitch-angle deflection of particles with mass m 
and charge ±e in fully-ionised hydrogen plasma (see ^.g|Emslie| 
[T9781 ). Thus, 


dv 


2n 


= -a- 


with 


a = 


(47rsoy 


rA 


-f nir 


meJnp 


where A is the Coulomb logarithm, and quasi-random changes 
of particle pitch-angles with probabilities, and, hence. 


v(t) -f Av = 


( 20 ) 


Corresponding quasi-random angular deflections are calculated 


as 


A6 = 1.3 ^/fiN VAt, 


( 21 ) 


where TV is a random number from normal distribution with 
the width of 1.0. Tests show that, statistically, this formula pro¬ 
vides a good approximation for pitch-angle diffusion as per 
I = // = cos6» = V||/v is the pitch- 


angle cosine ( [Landau 
written using Eqs.|20 


21 


1937). Therefore, collisional terms can be 


djd 

dt 

dt 


coll 


coll 


- V 




Av 

At 


dfi 

dt 


( 22 ) 

(23) 


- coll 


Initially, the box is randomly filled with 2^^ test-electrons 
which are uniformly distributed in space and in respect to pitch- 
angle cosine /i, while in respect to absolute velocity v particles 
have a Maxwellian distribution corresponding to T = 0.83 MK 
(the temperature of the coronal background; see Eig.[T]). Thermal 
bath conditions are applied on all six boundaries of the domain: 
particles are free to leave the domain, and for every particle leav¬ 
ing the domain one new particle is generated randomly from a 
Maxwellian distribution (at the thermal bath’s temperature) and 
is injected at the same point of the boundary. 

Hard X-ray bremstrahlu ng intensities are c alculated using 
the Kramers formula (see g.g jKontar et al.|20Q2| ): 




oo 

= e./ 


^ph 




where I(€ph) is the rate of emission of photons with energy €ph, e 
is an electron energy, fpeam is the differential beam electron den¬ 
sity (dnIdE), n is the space and time-varying electronic density 
of the ambient plasma, and Qo is the constant from the Kramers 
bremsstrahlung cross-section formula. 

The test particle approach is widely used to study particle 
acceleration in flares, but is not a fully self-consistent approach 
since it neglects the electromagnetic fields generated by the 
test particles themselves. However, in the configurations stud¬ 
ied here, the fraction of the released magnetic energy carried by 


the accelerated particles is relatively small - typically around 5% 
( jGordovskyy et al.||2014| ). Eurthermore, due to the fragmented 
nature of the current sheet, energetic electrons are farily uni¬ 
formly distributed within the loop volume, so we expect that 
their effects on the fields is fairly small. Hence, the test parti¬ 
cle approach is quite well justified. 


3. Magnetic fields, currents and hot plasma 

The kink instability in the present MHD simulations develops 
when the total twist reaches around 7 - 97i. It leads to the forma¬ 
tion of a helical current sheet around the twisted flux-rope (see 


Browning et al.j ( 

2008]); Hood et al. (2009|); Gordovskyy & 

Browning (|2012j); 

Gordovskyy et al. (20 14|);[Pinto et al.[([2015|). 


This, in turn, switches on the anomalous resistivity (where and 
when the current density becomes stronger than jcriu see Eq.[^, 
and leads to the start of the reconnection in the corona. There 
also are very strong currents near the loop foot-points, but the 
critical current threshold jcru in the lower atmosphere is much 
higher, and hence the resistivity there remains equal to the much 
lower background resistivity. Hence, the reconnection and cur¬ 
rent dissipation occur predominantly in the coronal part of a 
loop. The current density distribution evolves gradually into a 
more fragmented pattern thereafter, until it fades away. Magnetic 
reconnection occurs both within the twisted loop, and between 
twisted flux-rope field-lines and the ambient held. 

The reconnection episode is followed by a magnetic relax¬ 
ation phase lasting about 100 - 200 s in the smaller models, and 
about 300 - 400 s in the large-scale models. During this time, 
the magnetic energy is reduced by ~ 50 - 70 % of the free mag¬ 
netic energy in the system (Eig.|^. The amount of free magnetic 
energy is defined here as the difference between the magnetic en- 
ergy just before kink instability and the initial magnetic energy 
of the potential (untwisted) magnetic configuration. Most of this 
energy (around 80 - 90 %) is converted directly to thermal en¬ 
ergy due to Ohmic dissipation, while a small fraction is released 
in form of plasma kinetic energy. The latter is then also converted 
to thermal energy due to viscous effects, mostly around shocks. 
A more detailed description of the velocity held in reconnecting 
twisted loops is given by Gordovskyy, et al {submitted). 

Ohmic dissipation is much enhanced in the strong helical 
current, leading to plasma heating. Eigures and [7] show the 
spatial distribution of the plasma temperature and density at dif¬ 
ferent stages of the evolution of the system. Initially, the plasma 
temperature and density are everywhere equal to those of the 
initial gravitationally-stratified background {cf. Eig. [T]). During 
the thermal relaxation stage, thermal conduction increases the 
temperature and decreases the density in the transition region, 
while reducing the temperature and increasing the density in the 
lower corona. Eoot-point twisting leads to an increase in mag¬ 
netic pressure inside the twisted loop, making it expand. This 
expansion, in turn, results in a slight drop in plasma density and 
temperature inside the loops (by adiabatic expansion). Just be¬ 
fore the kink instability occurs, and near the loop-top, the tem¬ 
perature is about 40 - 50% lower than the temperature in the 
initial gravitationally-stratified atmosphere (see Eigures and 
[^a)), while the density is 50 - 60% lower (Eig.[7ja)). The pic¬ 
ture is more complicated closer to the foot-points: plasma around 
the expanding twisted flux-tube is being pushed vertically rather 
than horizontally due to the strong vertical density gradient. This 
results in accumulation of mass in a cylindrical shell around the 
loop’s legs. The density in these shells is about 2.5 times higher 
than outside, while the temperature is about 1.5-2 times higher 
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Fig. 8. Temporal evolution of the DEM for cases S, V and Y (left, middle and right columns, respectively). Time runs from the top to the bottom 
rows (as indicated in the inset legend on each individual panel, where time is indicated in units of to). For each one of the cases represented, the first 
row corresponds to the linear phase of the kink-instability, the second row represents the period of maximum thermal emission, and the third row 
the beginning of the relaxation phase. The upper tail which forms in the DEM distribution corresponds to the plasma heated after the development 
of the kink-instability in the magnetic flux-rope. This corresponds to a hot plasma component, which is particularly distinguishable in case V. The 
maximum temperature reached depends mostly on flux-rope’s magnetic field amplitude, approximately as ~ {cf. Pinto et al.|2015| . 


due to adiabatic compression. After the reconnection starts, the 
high temperature distributions near the loop-top resemble the 
current density distributions (e.g. |Hood et al.|20()9] ), starting with 
a helical shell, which then gradually becomes more fragmented 
and fades away. In contrast, closer to foot-points the hot plasma 
can be seen on the central axis of the loop. The density near the 
loop top evolves from a regular structure towards a fragmented 
one. Also, a substantial increase in plasma density after onset of 
reconnection can be observed in the loop’s legs, which is due to 
plasma evaporation from the lower atmospheric layers. 

Figureshows the DEM {T) of the coronal plasma at differ¬ 
ent instants of the simulation for three different models (S, V and 
Y). The DEM distribution extends into higher temperatures dur¬ 
ing the linear phase of the kink instability in all the models, form¬ 
ing a distinctive upper tail. The total emission measure of this 
hotter plasma component keeps increasing until the peak phase 
is reached, and a clear hot peak is sometimes visible (especially 
in model V). The bulk of this hotter plasma reaches temperatures 
about twice the temperature of the background plasma in these 
simulations. The DEM distribution flattens out afterwards, dur¬ 
ing the magnetic relaxation phase, but the maximum temperature 
remains higher than initially for a long period of time. Overall, 
the temporal evolution of the DEM analysed here is very similar 
to those discussed by |Pinto et al.| ( |2Ql3] ) in respect to a simpler 
model of a kink-unstable flaring loop. They related the forma¬ 
tion and the dissipation of hot flare plasma components in the 
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DEM to obse rvati ons by, e.g, |Reale et al.| ( |2009| ), [Battaglia"& 


|Kontar|p012| ) and jSylwester et al.| ( |2014| ). The only remarkable 
differences are that in our case the flare plasma reaches lower 
maximum temperatures and that the hot plasma component that 
forms during the course of the flare is less pronounced (because 
the ratio background to flare plasma is also higher here). 


4. SXR and HXR emission 

The combination of the variations in plasma density and tem¬ 
perature give rise to the soft X-ray thermal emission patterns 
displayed in Figs. The figure shows volume renderings of 
the thermal continuum emissivity at 2 keV at several instants 
{t = 34, 80, 136 and 227 s, respectively, from top to bottom) 
together with renderings of the plasma temperature. The left and 
right columns of each subplot show the same loop from differ¬ 
ent points of view (side/perspective and top views). The orange 
and blue lines represent magnetic field-lines rooted at the left 
and at the right foot-points, respectively. On the left subplot, 
the yellow/oran^ volumes represent the plasma emissivity e at 
2 keV (see Eq. [^, ranging from 6 = 5 x 10^^ (light yellow) to 
6 = 5 X 10^"^ erg • s“^ • cm“^ • Hz“^ (dark orange). On the right 
subplot, the light yellow/dark orange volumes represent temper¬ 
atures ranging from 1 MK to 2 MK. It is notorious that the con¬ 
tinuum thermal plasma emissivity does not exactly follow the 
evolution of the plasma temperature. The difference is due to the 









































Pinto et al.: Thermal & non-thermal emission from twisted coronal loops 


Emissivity (2 keV) 


Temperature 



5.0e+12 2.5e+14 5.0e+14 


Fig. 9. Volume renderings of the thermal continuum emissivity at 2 keV (left) and of the plasma temperature (right) at different instants for model 
S. The left and right columns of each plot show, respectively, a perspective view and a top view of the same loop. The colour-scales represent 
emissivities (e; see eq.|^ between 5 x 10^^ and 5 x 10^"^ erg • s"^ • cm"^ • Hz"^ and plasma temperatures between 1 and 2 MK. The orange and blue 
lines represent magnetic field-lines rooted, respectively, at the left and at the right foot-points. The instants represented are, from top to bottom, 
t - 34, 80, 136 and 227 s. The grey line delimits the bottom face of the numerical domain (a square of dimensions 20 x 20 Mm). 


variations in plasma density which result from the turbulent dy¬ 
namics of the kinking loop - note the dependence on density in 
Eq.0-, as was also the case in |Pinto et al.|p015] ). 


The first traces of thermal emission appear at mid-height and 
around the loop legs. These occur in the form of two separate 
sheet-like emission structures placed around the boundaries of 
the loop which are pushed against the background medium as 
the flux-ropes writhes. Thermal emission then fllls-up the whole 
flux-rope’s volume for a short period of time (at about the max¬ 
imum of emission). The total thermal emission fades out after¬ 
wards as the loop’s plasma cools down globally during the re¬ 
laxation phase. At this point, the thermal emission concentrates 
into a few separate field-aligned threads, most notably at mid¬ 
heights and crossing the top part of the loop (see, especially, the 
top views in Fig. |^. As the figure shows, thermal emission is 
only clearly discernible well past the onset of the kink insta¬ 
bility, when the flux-rope has already lost an important fraction 
of its initial twist. Therefore, the amount of twist which can be 
measured by tracing out the emission threads is much lower than 


the value of the initial magnetic twist. Both quantities differ by a 
factor 4 in our simulations (the emission threads show a twist of 
about In while the initial magnetic twist is ~ S/r). This situation 
is analogous to that described by [Pinto et al.| ( |2015] ) for a simpli¬ 
fied model of a kink unstable coronal loop (without large-scale 
curvature). Note that the thermal emission shows a structure with 
some writhe, which is in the sense of the original field-line twist. 
Fig, represents the thermal photon flux at the same photon 
energy (2 keV) and at the same times as Fig. but integrated 
along the line-of-sight for three different viewpoints (side, front 
and top of the loop), with the darker tones representing stronger 
photon flux. The contour lines show the locations of the peaks 
of non-thermal hard X-ray bremsstrahlung (HXR) emission at 
5 and 12 keV (blue and green contours, respectively, at 90% 
and 75% of the peak emission). The hard X-ray emission spa¬ 
tial profiles are smoothed with a Gaussian Alter with ~ 1 Mm 
width, roughly corresponding to the resolution of modern instru¬ 
ments (we are not trying to reproduce a particular instrument’s 
performance, however). It is interesting to note that, depending 
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on the line-of-sight (l.o.s), the thermal emission may appear as 
coming from an unique coronal loop or from two independent 
sources. This happens because the emission structures which are 
aligned with the l.o.s. are much enhanced in respect to the emis¬ 
sion structures oriented orthogonally, as the former occupies a 
shorter depth along the l.o.s than the latter. The top views (right¬ 
most columns) on Fig. and on Fig. [T^ show this effect very 
clearly. The emission threads which cross the loop above its apex 
are well visible in the volume renderings in Fig. [9 but nearly 
disappear when integrated along a vertical line of sight (Fig.p^. 
This effect further hinders the possibility of measuring the actual 
levels of magnetic field-line twist directly from the morphology 
of the thermal emission threads. 

Generally, the non-thermal electron behaviour we observe 
here is similar to that described in |Gordovskyy et al.| ( |2013| ). 
However, electrons appear to lose their energies faster. Most 
likely this happens because of the higher density closer to the 
footpoints. Thus, in the present model thermal conduction in the 
initial atmosphere leads to an increase of the density by factor 
of 2 around z ^ 3 - 5 Mm over footpoints and by factor of ~ 3 
around z 2 - 5 Mm in horizontal field (see Fig. [^. In addi¬ 
tion, dense plasma is also accumulated around loop’ legs dur¬ 
ing the twisting (Fig. [7]). Similarly, hard X-ray bremsstrahlung 
emission evolves slightly faster than in previous simulations. Al¬ 
though there are plenty of energetic electrons in all the coronal 
loop, the bulk of the emission is coming from very flat footpoint 
sources with the height of 1 - 2 Mm (comparable to the resolu¬ 
tion of these intensity maps). But the situation could, in princi¬ 
ple, be different in coronal loop configurations with higher coro¬ 
nal densities (HXR emission peaks could eventually appear at 
higher coronal heights, as is sometimes observed in solar flares). 
There HXR footpoint emission occasionally shows slight asym¬ 
metries, but these disappear quite quickly. On the overall, the 
HXR emission is essentially similar and simultaneous on both 
footpoints. The area of HXR footpoints initially increases with 
time due to the twisted loop reco nnection with ambient field (see 
[Gordovskyy & Browning|20T2| ). This increase is not monotonic 
and, at some point, the footpoint areas begin to shrink. This hap¬ 
pens due to two different factors: numerical undersampling due 
to the limited test-particle number, and modulation of the HXR 
intensities. Figure m show the spatial distribution of electrons 
with energies above 5 and 12 keV at different moments of the 
temporal evolution of model S (the same as in Fig.[T^. Electron 
distribution functions are calculated by reducing test-particles on 
a 40 X 40 X 40 X 40 [x, y, z, e] grid with uniform spatial step and 
exponential energy step to reduce undersampling at higher ener¬ 
gies. It can be seen that it is rather similar to the distribution of 
the parallel electric field in the domain (see Fig. [^, and hence 
that energetic electrons are distributed within most of the flar¬ 
ing loop volume, in spite of the HXR emission in Fig, being 
much more concentrated in the lower layers. Note that the HXR 
emission spatial patterns are not a direct reflection of the elec¬ 
tron acceleration mechanisms, but rather of the interaction be¬ 
tween the accelerated electrons with the highly stratified coronal 
plasma. A substantial amount of energetic particles accumulates, 
nonetheless, near the loop footpoints. This is because the mag¬ 
netic field convergence reduces the parallel component of the 
particle velocities. The magnetic trapping is even more effective 
due to collisional pitch-angle scattering moving particles out of 
the loss-cone. Since all the magnetic field lines of the twisted 
loop are tied to the lower boundary (i.e., the photosphere), virtu¬ 
ally all energetic particles leaving the domain do so through the 
lower boundary. 


Figure 12 shows the light-curves of the total thermal con¬ 
tinuum emission integrated in three photon energy bands be¬ 
tween 1 and 12 keV (see inset legend) and of the hard X-ray 
emission at 5 and 12 keV for models C, S and V. The light- 
curves are all normalised to their peak value for easier repre¬ 
sentation (absolute photon fluxes decrease with increasing pho¬ 
ton energy). The emission in lower energy bands is sustained 
for longer times (compared to non-thermal emission), as com¬ 
monly observed in flares. Photon fluxes at lower energy bands 
(e.g, 1-3 keV) systematically have a slow rise and decay, while 
flux at higher energy bands (e.g, 6 - 12 keV) always displays a 
more rapid evolution (as in |Pinto et al.|2015] ). The peak photon 
fluxes do not occur simultaneously at all photon energy bands: 
the peaks in higher energy light-curves tends to precede those on 
lower energy light-curves. This is particularly noticeable in mod¬ 
els S and V. A few post-flare oscillations are visible in all light- 
curves represented (hence, in an energy range between 1 and 
25 keV). As one would expect, non-thermal emission evolves 
much faster than thermal X-ray emission. In model S, for ex¬ 
ample, the HXR flux peaks at about r = 60 s, already after the 
saturation of the kink instability and the triggering of the fast 
reconnection, when the initial large-scale helical current sheet 
has already started fragmenting. This is visible by comparison 
with Figure the HXR peak occurs just in between the second 
and the third panels (t = 30 and 120 s, respectively). The cur¬ 
rent density distribution is indeed very much fragmented during 
this whole time interval, even though its large-scale spatial or¬ 
ganisation is a remnant of the initial helical current sheet. Both 5 
and 12 keV bremsstrahlung emission decay thereafter and nearly 
disappear after ~ 100 - 200 s, when thermal emission is still 
reaching its peak values. 

5. Discussion 

We study here the consequences of the triggering of the kink 
instability in twisted coronal loops in terms of X-ray emission 
(thermal continuum SXR and bremsstrahlung HXR) by means of 
numerical MHD simulations combined with test-particle meth¬ 
ods. The MHD model consists of a semi-circular kink-unstable 
coronal loop embedded in a stratified atmosphere composed of 
coronal and chromospheric layers. The resulting plasma den¬ 
sity and temperature are used to calculate emission measures 
and thermal continuum emission (soft X-rays), and the electric 
and magnetic fields are used to compute test-particle trajectories. 
Non-thermal (hard X-ray) emission is estimated based on the 
interaction between these test-particles and the stratified back¬ 
ground. Heating caused by these collisions and currents gen¬ 
erated by the accelerated particles are not fed back into the 
MHD simulation. That would require a different type of ap¬ 
proach based on hybrid methods, which is clearly beyond the 
scope of this manuscript. 

The vertical grid resolution-0.16 Mm and ~ 0.04 Mm 

for the "large" and "small" models, respectively - under-resolves 
the steep transition region (TR) between the chromosphere and 
the corona, even though they are at the forefront of current day 
simulations of the solar corona. For comparison, our vertical res¬ 
olution is similar or higher than the landmark three-dimensional 
sim ulation s by |Gudiksen & Nordlund| ( [2005] ) and |Bingert & Pe- 
ter ( 2Q11| ), respectively 0.15 Mm and 0.23 Mm, and approach the 
vertical resolution requirement at the TR of 20 km suggested by 
Klimchuk| ( [2006| ). Under-resolving the TR can have an impact 
on the properties of the simulated corona (density, temperature 
and/or heat fluxes), hence possibly affecting both the synthetic 
SXR and HXR emissions. In fact, [Bradshaw & Cargin] ( |2013| ) 
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Fig. 10. Thermal continuum photon flux emitted at 2 keV (colour-scale) and hard X-ray emission at 5 and 12 keV (blue and green contours, 
respectively, at 90% and 75% of the peak emission at those photon energies) integrated along the line-of-sight in three different directions (from 
left to right: side, edge and top views), for our reference case (model S). The thermal and non-thermal photon fluxes are normalised to arbitrary 
units for easier comparison. The axis are in units of Mm. The instants represented are, from top to bottom, t - 34, 80, 136 and 227 s (as in Fig. 
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have shown that under-resolving the transition region may sig¬ 
nificantly underestimate the peak density of coronal loops under¬ 
going impulsive heating episodes (and being loaded with plasma 
by means of chromospheric evaporation), but without much ef¬ 
fect on the plasma temperatures they reach. In spite of these limi¬ 
tations, our results are consistent with the main common features 


of observed solar flares, and furthermore allow us to clearly es¬ 
tablish links between different types of X-ray emission. 

Our models represent confined (non-eruptive) flares in single 
loops. The cases studied here represent smaller - and also thus 
more frequently-occurring - events. Larger energy releases can 
be obtained by considering stronger magnetic fields and larger 
loops (Gordovskyy et al 2015 submitted). However, whilst the 
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a) 




b) c) 


Fig. 11. Spatial distribution of energetic electrons in the model S seen from the loop side (cf. the first column of Fig.[^ at different instants. 
Left column (a), middle (b) and right (c) columns correspond to 35 s, 81 s, 104 s after onset of the kink instability. The axis indicate distances in 
Mm. Red and green colours indicate high and low electron densities, respectively. Upper and lower panels correspond to particles with energies 
> 5 keV and >12 keV, respectively. 


kink instability in a loop may play a role in the largest flares 
(e.g. by triggering a filament eruption), such events would in¬ 
volve larger-scale and more complex field configurations than 
considered here. Our results provide a prediction of what is ex¬ 
pected to be observed in kink-unstable loops, and thus new ob¬ 
servational campaigns are required to look for such signatures. 
Detailed comparison between simulations and observations is a 
task for future work. However, much work already published 
provides evidence of events which appear to be consistent with 
our models. The simulated loops rise somewhat during the main 
reconnection phase, then contract back down, so may be ob¬ 
served as "failed eruptions", which appear sometimes to be as¬ 
sociated with kink instability (e.g. | Alexander et al.||2006HSong| 
let al.|2014 ). Flares can be directly associated with apparent in¬ 
stability of observed twisted flux ropes (g.g jWang et al.||2015| ). 
|Liu et al.] ( |2013| ) analyse a flare within a loop, suggesting that 
the x-ray emission is generated by fragmented current sheets in 
the coronal part of the loop. Furthermore, there is evidence for 
an association of a p eriod of helicity build up prior to flare onset 
( |Reinard et al.|201()| ), consistent with our scenario. 

Our results suggest that highly twisted coronal loops may 
be more abundant than what current observations seem to indi¬ 
cate. In fact, individual magnetic fleld-lines in flaring coronal 
loops will only be visible as emission threads if (and only if) 
the enclosed plasma has already been considerably heated in re¬ 
spect to the background plasma during the course of the flare. At 
that point, magnetic fleld-lines twisted well above the threshold 
for the kink instability will already have lost an important frac¬ 
tion of their initial twist. The effect is here even stronger than 


in the simulations reported by |Pinto et ^ ( |2015| ), in that part 
of the emission threads nearly disappear in the line-of-sight pro¬ 
jections. Comparing Figure which shows the photon fluxes 
integrated along the line-of-sight for the three-dimensional emis¬ 
sion structures shown in Fig.|^ makes this particularly evident. 
The threads which cross the loop above its apex nearly disappear 
when integrated along a vertical line of sight, for example. These 
emission patterns could even be interpreted as coming from a 
coronal structure with a slight large-scale torsion but with no 
twist at ah. This result removes one of the main obstacles faced 
by the kink instability scenario for the theoretical understanding 
of the triggering and evolution of solar flares. 

The temporal ordering of the soft and hard X-ray light- 
curves we obtained are consistent with those typically observed. 
The light-curves in Fig. show similarities with GOES and 
RHESSI light-curves, such as those reported by jSylwester et al. 
( |2014| ) for an Ml.O class flare. Lower energy thermal (SXR) 
light-curves grow and decay more slowly than higher energy 
thermal light-curves ( |Pinto et ST]|2015| in a way similar to that 
in). Eurthermore, the hard X-ray (HXR) light-curves due to the 
bremsstrahlung of non-thermal electrons reach their peak value 
when the growth-rate of the SXR curves are maximal. This be¬ 
haviour is consistent across ah the models we tested (C, S, V 
and Y), and approximately reproduces the so-called Neupert ef- 
fect.The most commonly accepted interpretation of this effect 
invokes heating of chromospheric plasma by the supra-thermal 
electron beams which produce HXR emission, followed by the 
up-flow of hot plasma into the coronal loops and subsequent 
SXR emission. However, our simulations suggest a different 
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Fig. 12. The top row shows light-curves of the soft X-ray emission (black and red lines) and hard X-ray emission (blue and green lines) in models 
C, S and V at different photon energy bands (see the inset legend). The bottom row shows the temporal derivative of the soft X-ray lightcurves 
together with the hard X-ray lightcurves during the first 150 s for the same models. The curves are all normalised to their peak value for easier 
representation. The peaks in the higher energy SXR light-curves (red lines) tend to precede the lower energy ones (black lines); this especially 
visible in models C and S. The decay is also faster at higher SXR photon energy bands. The HXR light-curves (blue and green lines) tend to peak 
at the moment when the SXR growth rate is at its highest (this is especially visible in model S), such that the temporal derivatives of the SXR 
curves correlate well with the corresponding HXR curves. A few post-flare oscillations are visible in all the SXR light-curves represented 


cause for this effect, as the electron beams we consider are not 
allowed to heat up the background plasma by construction (the 
supra-thermal electrons are test-particles). Here, the instant at 
which a given HXR light-curve peaks is very close to the time 
when the current sheet that forms around the twisted coronal 
loop reaches its maximum amplitude and spatial extent (i.e, just 
before the main magnetic reconnection episode stars). Particle 
acceleration is maximal at this instant, and so are the flux of 
particle reaching the dense lower layers of the atmosphere and 
the consequent HXR photon flux. Then, actual plasma heating 
leading to SXR emission keeps occurring afterwards during the 
relaxation phase (even though the ohmic heating is proportional 
to 777 ^, and hence maximal at that instant). This result is in agree¬ 
ment with some of the conclusions of |Veronig et al.| ( [2005| ), 
namely that fast electron beams may not be the main source of 
SXR plasma supply and heating, and that other heating sources 
are necessary to explain Neupert effect quantitatively. 

The ensemble of results presented in this manuscript sup¬ 
port the view that plasma heating and particle acceleration pro¬ 
cesses in solar flares may occur in large volumes of the respec¬ 
tive flaring loops rather than in a small volumes located near their 
apexes. The flare scenario based on the triggering of the kink- 
instability of coronal loops with magnetic helicity (twist) natu¬ 
rally provides the physical mechanisms necessary for volume- 
distributed heating and acceleration. 


6. Summary 

The aim of this study is to determine observational manifesta¬ 
tions of solar flares occurring in kink-unstable twisted coronal 
loops. Specifically, we focus on thermal soft X-ray (SXR) con¬ 
tinuum emission together with non-thermal hard X-ray (HXR) 
bremsstrahlung emission deduced from a set of MHD simula¬ 
tions. The thermal SXR emission is calculated from the distribu¬ 
tions of density and temperature of the compressible plasma, and 
the HXR emission is derived from the interaction between test- 
particles and the time-evolving magneto-hydrodynamical strati¬ 
fied background. 

Our results show that the geometry of the emission produced 
during flares does not necessarily correspond to that of the un¬ 
derlying magnetic field. In particular, the twist angles revealed 
by the emission threads (soft X-ray thermal emission; SXR) 
are consistently lower than the actual field-line twist present 
at the onset of the kink-instability. Hard X-ray (HXR) emis¬ 
sion due to the collisions of energetic electrons with the strat¬ 
ified background are concentrated at the loop foot-points (quasi- 
simultaneously), event though the electrons are accelerated ev¬ 
erywhere within the coronal volume of the loop. The maximum 
of HXR emission consistently precedes that of SXR emission 
The HXR light-curve being approximately proportional to the 
temporal derivative of the SXR light-curve, hence reproducing 
the so-called Neupert effect. 
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Future work will cover other aspects of thermal emission in 
flares such as thermal emission in the extreme ultraviolet en¬ 
ergy bands (continuum and line emission), and also properties 
of micro-wave gyro-synchrotron emission in flares triggered the 
kink instability. 
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